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Abstract 

We calculate the partition functions of QCD in two dimensions on a cylin- 
der and on a torus in the gauge OqAq = by integrating explicitly over the 
non zero modes of the Fourier expansion in the periodic time variable. The 
result is a one dimensional Kazakov-Migdal matrix model with eigenvalues 
on a circle rather than on a line. We prove that our result coincides with the 
standard expansion in representations of the gauge group. This involves a 
non trivial modular transformation from an expansion in exponentials of to 
one in exponentials of l/g'^- Finally we argue that the states of the U{N) or 
SU{N) partition function can be interpreted as a gas of N free fermions, and 
the grand canonical partition function of such ensemble is given explicitly. 
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1 Introduction 



The last year has seen a revival of attempts to understand non-perturbative QCD 
using the 1/N expansion . First, Kazakov and Migdal (KM) proposed a model 
with an adjoint multiplet of scalar fields coupled to external gauge fields, which 
turned out |0| to be exactly solvable in the large N limit. The model was originally 
proposed as a way of finding the master field that would solve four dimensional QCD 
in the large N limit. It is not clear yet whether this goal will be achieved or not, as 
it is almost certain that some modifications to the original model are needed for it 
to admit the QCD fixed point in the continuum limit, but the proposal has already 
generated a considerable literature (see for example and references therein ). 

It has been pointed out that the ci-dimensional KM model with a quadratic 
potential describes the high temperature limit of pure QCD m d + 1 dimensions; 
however, for d> 1 such theory does not have a continuum limit unless higher order 
terms are added to the potential. On one hand, this is in agreement with known 
results in dimensional reduction of finite temperature QCD, but, on the other hand, 
higher terms in the potential are difficult to control and their arbitrariness limits 
the predictive power of the model. 

Whatever its relationship with four- dimensional QCD, it is remarkable that the 
KM model is the only known example of a matrix model, hence of a string theory, 
that is solvable in more than one dimension in the large limit. The question of 
what kind of string theory corresponds to the KM model has also been addressed 
recently 0. It turns out that the KM model describes, at least in the large N limit, 
a string theory with infinite string tension, where the string has collapsed into a 
branched polymer. This might explain the solvability of the KM model even in 
dimensions higher than one. 

More recently, in quite a different context. Gross and Taylor proved that 
QCD in two dimensions (QCD2) can be interpreted as a string theory in the large 
N limit, by showing that the coefficients of the expansion of the partition function 
of QCD2 in power series of can be interpreted in terms of mappings from a two 
dimensional surface onto a two dimensional target space. A lot of information on 
the underlying string theory is encoded in the coefficients of this expansion, but the 
lack of a prescription for the string action and for the calculation of the string path 
integral limits the efforts to gather non-perturbative information or extrapolate to 
higher dimensions. 

In this paper, we proceed in a different direction in the understanding of the 
string theory underlying QCD2, by proving that QCD2 on a cylinder and on a torus 
is described by a matrix model which is exactly a one dimensional KM model on 
a circle with one extra condition: that each eigenvalue of the scalar field in the 
KM model be defined modulo integer multiples of 2tt. In short, QCD2 on these 
surfaces is a periodic one- dimensional KM model, and dimensional reduction (or 
equivalently the high temperature expansion) is exact in this case. 

We emphasize that we are not taking the large limit, and our proof is valid 
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for any N. Furthermore, the method we employ leads us to find expressions for the 
partition functions that involve exponentials of the inverse gauge coupling, unlike 
the usual character expansions 0, ||, 0- The existence of such expressions is a con- 
sequence of the well-known fact that the partition functions of non-abelian gauge 
theories on two-dimensional surfaces with boundaries are kernels of the heat equa- 
tion on the gauge group manifold, as expected from the path integral formulation. 
Such kernels, at least for U{N) and SU{N), are expected to admit representations 
in terms of periodic gaussians in the invariant angles (sometimes called Weyl angu- 
lar parameters, they are just i times the logarithm of the eigenvalues of the unitary 
matrix in the fundamental representation). Indeed, such a representation has been 



known for a while in the case of the partition function on the disk ||10| , \L1 . 

By a careful gauge-fixing procedure, we are able to find a similar representation 
for the kernel associated with the cylinder, and this leads to the partition function on 
the torus via identification of the boundaries and group integration. The invariant 
angles that parameterize our representation are precisely the eigenvalues of the 
KM model, and are directly related to the Polyakov loops around the cylinder, as 
expected from [Q. 

After deriving the partition functions, we are able to explicitly prove the equiva- 
lence of our expressions with the corresponding character expansions, by performing 
a modular transformation that generalizes to non-abelian gauge theories the mod- 
ular inversion of Jacobi's ^3, which in fact is the partition function of QED on a 
torus. 

As a result of our analysis we also find a generating function for the partition 
functions of U{N) [SU{N)] . As expected from the identification of QCD2 with a 
one dimensional matrix model, this grand canonical partition function is given in 
terms of infinite products of free fermion partition functions. 

The paper is organized as follows: in section 2, after a few general remarks, 
we calculate the partition function on the cylinder with an appropriate choice of 
gauge fixing, and we get a representation in terms of angles that coincides with 
the announced KM model. In section 3 we proceed to reconstruct the character 
expansion on the cylinder by a modular transformation. In section 4 we calculate the 
partition function on the torus and we find the generating function that reproduces 
the partition functions of SU{N) {U{N)) for any N. In section 5 we present our 
conclusions, and finally in the appendix we present an alternative derivation of the 
kernel on the cylinder, following closely the techniques applied to hermitian matrix 
models in . 



2 The Partition Function on the Cyhnder 

It is well known by now that two dimensional QCD defined on a manifold A4 of 
genus G and with a metric g^i, is exactly solvable. The partition function is given 
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by 



R 



where the sum is over all equivalence classes of irreducible representations i?, dR 
is their dimension and C2{R) is the quadratic Casimir in the representation R. 
Similarly, the heat kernel defined by a surface of genus G and n boundaries is given 

by ii 

JCcAai, ...,gn,N,A) = Y: dl'^^-^XRigi) ■ ■ ■ Xi?(^?n)e-^-^^'^^(^\ (2) 

R 

where gi are the Wilson loops evaluated along the boundaries, and xr denotes the 
Weyl character of the representation R. For dimensional reasons, and because of 
the invariance of the action under area preserving diffeomorphisms, eqs. (|l|) and 

depend only on the variable g^A. We will henceforth denote this variable by t. 

The heat kernels in eq. (§) are class functions, and therefore they must admit 
a representation in terms of the eigenvalues of the group elements g^ in the funda- 
mental representation, which are pure phases for SU{N) and U{N). An explicit 
expression in terms of these phases was found by Menotti and Onofri for the 
simplest case, the disk. It is of the form 

r,,|tloo *<iii 2 sm 2 [0i - 0, + 2'K{k - /,)] 



{/,}=-cxD i<i 

X exp 



1 ^ 
i=i 



(3) 



In order to formulate QCD2 as a KM model, we will construct a similar repre- 
sentation for the kernel defined by the cylinder (/Co,2); which we will denote simply 
by /C2. It will be convenient to work in first order formalism, where the action is 
given by 

S[N,t) = - r dxdr Tr[F2 - iFf{A)V (4) 
t Jo 

where 

fiA) = doA,-d,Ao-t[Ao,A,], (5) 

while F and A are considered as independent fields. By using the invariance under 
area preserving diffeomorphisms of the original action, we have restricted the metric 
in eq. to be fiat and rescaled space-time coordinates (x, r) so that they both 
range in the interval (0, 27r). The variable t appears then explicitly at the r.h.s. 
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In the case of a cylinder we have to fix the value of the Polyakov loop at the 
two boundaries (say at x = and x = 2tt) by introducing suitable delta functions : 



IC2{9i,92,t) = JvA^VFexpi^--TTj^\xdT[F'-iFf{A)]^x 
X 6 {W{0),g,) 6 {W{2n),g2) H9i)H92), 



(6) 



where 



and 6{g, h) denotes the conjugation invariant delta function on the group manifold, 
defined by 

5{g^h) = f dU6{UgU^h). (8) 



The factors ipioi) V'(5'2) are just normalization factors; they depend only on the 
eigenvalues of gi and g2 and they will be chosen so that the sewing of two cylinders 
corresponds to just a group integration: 



J ^2(5-1, g, )/C2(fi'^ g2)dg = /C2(fi'i, fi'2)- 



(9) 



All fields in eq. (|^) are periodic with period 271 in r and can be expanded in Fourier 
series as 



(10) 



All the Fourier modes with n ^ can be integrated away in eq. (^ with a suitable 
gauge choice, namely 





{x 


r) ] 


+00 


" Bn 


i^) 1 


A, 


{x 


r) 


= E 


An 




F{ 




r) 


n=— 00 


Fn 


» . 



doAoix,T)=0 ^ 5„(x) = 0, n^O. 
In this gauge the kernel becomes 

)C2{gi,g2,t) = fvBol[l^AnVFne-^^'C''^^(^'^'^MAFp 

X S iWiO),g,) S iW{27r), (72) 

where 

+00 

C= {FnF_n-nFnA_n + tSn,oFodiBo-Fn[Bo,A_n]}, 

n=— 00 

while Ai?p is the Faddeev- Popov determinant 

^6[nen - [Bo,en\] rs 



X 



App = Y[ det- 



6e 



n,tl 



and W{x) is now given by 



W{x) = e 



2inBo{x) 



(12) 
(13) 

(14) 
(15) 



5 



At this point the functional integration over An {n ^ 0) leads to a product of delta 
functions : 

n 5 {nFn - [i?o, Fn]) = -^X[ 5{Fn). (16) 



By inserting this equation in eq. ([1^) we obtain 

X 5 (ly (0), 5 (ly (277), (?2) ^(^7i)^(^72) 

X 5(iy(0),^7i)5(iy(27r),(72)^(^?i)^(^72). (17) 

where we have replaced i?o, aiid Fq with i?, A and F. This result can be in- 
terpreted as follows: the action is a KM model in one continuum dimension with 
boundary conditions depending on e^™^ rather than B\, this implies that if, for 
instance, -B(O) satisfies the boundary conditions imposed by the delta functions at 
X = 0, then the latter are satisfied by any other matrix whose set of eigenvalues 
coincides with the one of -B(O) modulo integer numbers. It is indeed possible to re- 
strict all the eigenvalues of Ao(x, r) = B{x) to be in the interval [0, 1) by a further 
gauge choice. In fact if Aq{x,t) = ?7(x)diag (Aj(a;)) W{x) then the gauge transfor- 
mation generated by the unitary matrix h{x, r) = ?7(x)diag (e*'^*'^) W{x) is periodic 
in r with period 27r, preserves the gauge choice (p!TD, and just amounts to a shift 
of all the eigenvalues Xi{x) by the integer numbers n^. It is clear though that to 
avoid discontinuities in the gauge transformation the integers rij in h{x,t) have to 
be X independent , so the eigenvalues Aj(x) can be restricted to the interval [0, 1) 
only for one particular value of x, say x = 0. By means of techniques which are by 
now standard in matrix models, the matrix B can be diagonalized and the action 
written in terms of its eigenvalues Xi{ 

JC2{gi,g2,t) = lvK{x)A{X{0))A{X{2n))expl^-^ l^\xdK{x)d\\x)^x 

X / dUidU2S [t/ie2^"^(°)f/|e-*^l 6 \U2e^'^^^^''^Ule-''^] V^(^)^(0),(18) 



where we inserted the explicit expression for the 5 functions, we denoted by A(0), 
A(27r), 6, (j) the diagonal matrices whose eigenvalues are respectively Ai(0), Ai(27r), 
9i, and 0,, and we chose e~*^ [resp e"""^] to denote the diagonal form of gi [resp of 
g2\. A(A) is the Vandermonde determinant of a hermitian matrix with eigenvalues 
\. 

The integrals over the unitary matrices Ui and U2 can be calculated from the 
identity 



exp 



■;5Tr(\/ + _ 2)1 = ^rf^/5-^V2^^(\/)(l + 0(l//3)) 



R 



r^'/^5(y), (19) 
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where is a unitary matrix. 

Inserting this expression of S{V) in eq. 
to calculate the integral for large /3f|, we find 



and using the saddle point method 



dU6 ( Ue'^U^e'^ 



+00 

E E 

P {n,}=- 



.l)P+iN-l)j:,n, 



N 



where P denotes a permutation of indices, 

J{a) 



n5((T,+ePO) + 27™,), (20) 



[]2sin 

i<j 



(21) 



is the Vandermonde determinant for a unitary matrix, and we made use of the fact 
that J{a + 27rn) = (-1)^'^-')^. "V(a). 

The quadratic functional integral over the eigenvalues can be performed by a 
^-function regularization of the divergences, leading to 



VXi^x) exp 

1 \^/2 



1 



/ dxY^{dK{x)y 

Jo ~1 



-V 



exp 



N 



l5:(A,(0)-A,(27r))^ 
■^'^ 1=1 



(22) 



The result at the r.h.s. of eq. ( p2D is the same that we would have obtained if 
we had started from a KM model on a one dimensional lattice with an arbitrary 
number of sites in the interval [0,27r). The result is independent of the number of 
sites due to the remarkable scaling properties of the KM model, already noticed in 
T^. By substituting eq. 



and eq. ( POD into eq 

^-N/2 



one finally obtains 



exp 



N 



i=l 



(23) 



Here the normalization factor has also been determined. It is given by il){g) = 
J{g) / A{g), where A{g) = ni<j(^i — ^j) if the eigenvalues of g are denoted by e*^\ 
For the gauge group U{N) the eigenvalues (pi, 6.1 and the integers U are un- 
constrained, whereas for SU{N) we can choose J2i 4>i = = 0, and this in turn 
constrains the integers k to obey J2i h = 0. In fact if we choose the gauge Aj(0) = 9i, 
as previously discussed, then the fact that Aj(27r) = (pi + 27r/j, together with the 
continuity of the function J2iK{x), implies J2ik = 0. Thus for SU{N) the sign 
factor in each term of eq. (p3|) is just given by (— l)"^- Note also that for SU{N) 
the factor t"^/^ must be replaced by 

""^ These generalizations of the Harish-Chandra integral are discussed in section 3 of ref. jl^ 
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3 Modular Inversion of the Kernel on the Cylin- 
der 



The expression for the kernel on a cyhnder obtained in the previous section, as 
announced, involves exponentials in unlike eq. As a preliminary check on 
eq. p3|), one can easily verify that by taking the limits ^ — or — >^ one obtains 
the kernel on the disk, eq. (||). 

To verify that the two expressions are related by a modular transformation, let 
us consider now eq. (|^) for the gauge group SU{N) : 



f2lT 

rrl JO 



1 t{l-A^)/2 
/ 1 ^ r 

exp {^P'{i) - ^P(i) + 27r/i 



p+p' 



2i(3tli 



(24) 



The lagrange multiplier (3 has been introduced to impose the condition I] = 0, 
and a redundant double sum over permutations has replaced the simple sum of 
eq. 

By completing the square, eq. (El) can be rewritten as 



lC2igi,g2,t) 



1 ^ 



27rN\ 



exp 



E 

p,p' 

271 



J{0)J{4>) jj^ 



E/ 



p+p' 



2 N 



X exp 



t 

'"8^ 



Eh 



(l>p'{i) — Op(i) . (3 



1=1 

2\ 



2n 



47r2 



X 



(25) 



so that one can use the well known modular transformation of the function 6- 



+00 



E 



l=—oo 



At 



+00 



exp (-47r^n^tj exp(27rm6')(47rt) 



a/2 



to obtain 



E 



^1/2 



2vriV! ^, J{e)J{<j>) ^ 



E / df3{-V 

{rii} 



p+p' 



N 



N 



X exp 

_ 

~ 2ttN\ 



E^? + 2vr.E 



0p'{j) — Op{i)_ _ . (3 ^ 



i=l 



1=1 



27r 



E 



47r2 



X 



Ntl3^ 
87r2 



3) 



(26) 
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X ~, — 7T _ , „. X 



J(0) 



2-K 



d[3 exp 



Nt I 



(27) 



/5 



27r 



E 



n,; 



The sum over {nj} can now be restricted to the region ni > 77-2 > • ■ ■ > riAr, at the 
expense of the factor 1/A^!. One can further notice that both the first exponential 
and the determinants are invariant if all the rij's are shifted by the same constant. 
One is then lead to define ri = rti — nj^ — N terms of which the kernel becomes 



^2(91, 92, t) 



2ti 



E exp 

ri>r'2>...>rjv=— 



X 



det{e"'"^j} det{e-*'"'^^} 



J(0) 

+00 

E 



/3i 



m 

dp exp 



X 



(2J 



Nt 
"8^ 



P 



27V 



E^ 



where Pq = —27rnN — 2tiN = Pi — 27r. The sum over simply reconstructs the 
gaussian integral from —00 to +00, that can be trivially calculated. It is easy to 
recognize now that in eq. (28) the exponential is related to the quadratic Casimir 
of the representation R whose Young tableaux has rows of length fi = ri + i — N. 
In fact 



E^ 



1 

iV 



E^O = C2(i?) + ^iv(iv2 - 1). 



(29) 



The last term can be interpreted as a zero point energy, namely the energy of 
the lowest representation. It is proportional to the scalar curvature of the group 
manifold, and to the modulus squared of the vector obtained by summing the 
positive roots of SU{N). On the other hand we have also 



det{e*'""^^} 



N(N-l) 
-t) 2 



(30) 



By substituting (^) and (|30D into eq. (28) we obtain the modular inversion for the 
SU{N) kernel on the cylinder, which reads 



exp l--N{N' - 1) E^^P -o^2(i?) XRhe)xR{<P) 



24 

It 

^1/2^ ^(l-iV)/2 ^ 

4 J ^J(^)J(0)|^"^P 



1 ^ 



P{i) 



27tL 



(31) 



where the constraints J2i h = J2i 



J2i = are understood. 
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4 The Partition Function On the Torus 



The kernel on a cylinder obtained in section 2 section allows one to calculate the 
partition function of QCD2 on a torus by simply sewing together the two ends of 
the cylinder, according to 

ZG=i{N,t) = J dgJC2{g,g,t) = J^'^ l[d4>,j\(P)]C2{^,^,t). (32) 

This was done in ref. by using eq. (H) with G = 0, n = 2 and the orthogonality 
properties of the characters. By repeating the same calculation with our expression 
(p^) for the kernel we are able to write Zg=i{N ,t) in terms of 9 functions, whose 
behaviour under modular transformations is well known. The result is particularly 
simple in the case of the group U (N) where in the expression of the kernel, eq. 
the integers /j and the angles 6i and (pi are unconstrained. 
In this case the partition function is given by 

^ X 



/ 1 Af \ 

X E(-l)^''"'^^^'^exp --^(0,-0p(,) + 27r/,f . (33) 

{h} V ^^i=l / 



By using the modular transformation eq. ( PB] ) for the theta function the expo- 
nentials can be made linear in the angles (pi : 



r2n N 



ZG=i{N,t) = [ ^ Yld<PiJ2i-^'^ 

X ^^P ( ~| (^i " + iJ2i'>^i- ^n) {(pi - <Pp{i))] , (34) 
{m} \ i=i i / 

where 6n is for odd and 1/2 for even N. The r.h.s. of eq. (153) can be computed 
by the same method employed in for the one-dimensional KM model, that is 
by decomposing each permutation into its cycles. Then 



/ N \ N 

hi,...,h]sf \r=l / r=l 



where is the multiplicity of a cycle of length r in a given permutation and the 
sum over all permutations is reproduced by summing over the 's with the correct 



^This partition function has a different zero point energy compared to the one defined for 
instance in ref. P] corresponding to the overaU factor exp[| ^'^^~^^ \ 
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combinatorial factors. Fr is the contribution from a cycle of length r and it is given 
by 



2tt 



I d(j)i . . . d(f)r exp 

{rii} 







^ i=l i=l 



= (2vrr E e-^("-^-)^ = (2vrr^.(^)(0,r = .|^), (36) 

n=-oo ^ ^"^ 

where 9„(^n) denotes Jacobi's 6*2 for even and 6*3 for N odd. 

This result can be expressed in a rather elegant and interesting form if one 
consider the grand-canonical partition function 

ZG=M = Y.ZG=i{N,t)q^ (37) 

TV 



The even and odd parts of ZG=i{q) have to be computed separately. From eq. ( P7| ) 
we have 

miiQ) = ^exp|f:^^2fo,r = z^)(2vrg)4 + (g<^-g), 



2 r ^ V ' 27r 

' 00 / -i\r 



2G=i{q) = ^exp|E^ll^e3(0,r = z^)(27rg)^|-(g^-g). (38) 

The sum over r in the exponents can be performed if one replaces the theta 
functions with their expressions as infinite sums, finally leading to the result 

^G=i(?) = -n l + 2vrge-^("+^) +(g--g), 

^ n=0 ^ ^ 

ZttM = -(l + 27rg)n(l + 27rge-|"') - (q ^ -q) . (39) 

^ n=l 

This partition function has some similarity with the one of the one-dimensional 



KM model |]T^, in the sense that they both describe a gas of free fermions. However, 
in this case there are some peculiarities. The integers n can be interpreted as the 
discretized momenta associated with the winding of the eigenvalues, therefore the 
energy levels are proportional to n"^, rather than n. For n 7^ two degenerate 
levels are present, corresponding to winding in opposite directions. On the other 
hand, the integers U in eq. (|33|) are the winding numbers, namely the discretized 
coordinates of such modes, and the modular transformation allows us to go from 
the coordinate to the momentum representation. Further, the energy levels depend 
on the parity of the total number of fermions, which is connected with the fact 
that a fermionic wave function picks up a different sign according to the number of 
other fermions encountered in the winding. This is the meaning of the phase factor 
in eq. 
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We also emphasize that Zc=i{q) has been defined keeping t = g'^A independent 
of N, so eq. ( |39|) cannot be used as such to calculate the large limit of Zg=i{N, t). 

At least for odd, eq. (|38D has a simple connection with QED2, as one can see 
by remembering that 63(0, i t/2n) is the partition function of QED2 on a surface of 
area A 0. Consider now QCD2 on a torus and go into the gauge where F{x,t) is 
diagonal. An easy calculation shows that the Vandermonde determinants resulting 
from diagonalizing F{x, r) are exactly cancelled by the functional integral over the 
off-diagonals elements of A^[x,t). This leaves one apparently with copies of 
QED. However, as one goes with continuity from x = to x = 27r the eigenvalues 
of F will be permuted in an arbitrary way, and the sum over all permutations is 
needed to recover the whole partition function. One should notice finally that, in 
the decomposition of each permutation into irreducible cycles, each cycle of length 
r corresponds to a QED2 where we have to wind r times in the x direction before 
we go back to the original point, so in fact a QED2 on a torus of area rA. 

Consider now the case of SU{N). The partition function for the group SU{N) 
is the same as eq. (|55|), but with A^ replaced by A^ — 1 and the constraint = 
J2h = 0- Thus 



1~N 



Air 



6 



V Jo 2vr^ 



X 



f 1 ^ 2^1 

exp |-— ^ (0i - + 27rZi)^ + ipYjiy 



(40) 



The sum over the U can be done by reconstructing the square and using eq. (p6D , 
in analogy to what was done in in the case of kernel of the cylinder. At this point 
the integral over the (^j is trivial and gives 



. ^ 1/2 N 

ii) EE(-i)'n^»...„ 

K} P i=i 



exp 



■ N 
i=l 



(E^» 

A^ 



X / d[3 exp 
/o 



Nt 



(41) 



It is not difficult now to show that this partition function coincides (apart from the 
zero point energy) with the one of ref . 0, ||, H , 



7C2(i?) 



(42) 



R 



In fact the determinant of 5 functions in eq. (^T]) just forbids any couple of n^'s to 
take the same value and the integral over f3 becomes trivial after the shift of the 
nj's leading to eqs. (27) and (28). 

Starting from eq. (|4l|) it is possible, following the steps leading to eqs. (|37D and 
(p9|), to write in compact form the grand canonical partition function for SU{N), 
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as 



dl^ n (^l + ge-^("-^^) j . (43) 



t /■27r 

47r. 

The interpretation of this formula is the following: the constraint Yl,ih = ^ means 
that for SU {N) the wave function of the center of mass of the eigenvalues is com- 
pletely localized, and therefore the corresponding momentum is undetermined. This 
is the reason of the integration over /3 in eq. p3|), and of the invariance of the 
Casimir in eq. ( PP] ) under common shifts of the integers labelling the representa- 
tions. 



5 Conclusions 

We have shown in this paper, by means of dimensional reduction techniques, that 
QCD2 on a cylinder or a torus is exactly a one dimensional matrix model of the 
type proposed by Kazakov and Migdal, with the substantial new feature that the 
eigenvalues of the matter fields live on a circle. This is in agreement with the 
interpretation of the scalar fields as the logarithm of the Polyakov loop in the 
compactified dimension. We prove also that the fundamental constituents of the 
theory are free fermion excitations corresponding to the winding of the eigenvalues 
around their target space, with energy growing quadratically with the winding 
number. As shown by the standard form of the partition functions, eq. the 
states of the theory are labelled by the irreducible representations of U{N) (or 
SU{N)), and are states of such free fermions. 

The interpretation of QCD2 as a matrix model is natural in a formulation where 
the kernel on the cylinder and the partition function on the torus are expressed as 
an expansion in exponentials of |^ rather than in exponentials of g'^A as in the 
representation expansion. From the mathematical point of view the bridge, between 
the two representations is a modular transformation that generalizes to the cylinder 



and to the torus the result already known for the disk |jTO| . 

It would be very interesting to generalize this result to kernels with a higher 
number of entries and to arbitrary genus, that is to find for 

JCcAgu 92,..., 9n, iV, t) = E dl-'''-''XR{9i) . . . Xii(^?n)e-^^^(^) (44) 

R 

a representation in terms of exponentials of ^j^. This would presumably amount 
to describe such kernels in terms of one dimensional matrix models whose target 
space includes both branching points and loops. There is one obvious difficulty in 
such program: by sewing two kernels together one should obtain a new kernel that 
depends only on the total area. This is quite natural in (0), as the exponentials are 
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additive in the area, but in principle much more difficult with expressions containing 



exponentials of 




6 Appendix A 

In this appendix we will provide an alternative derivation of the SU {N) partition 
function on the cylinder, which follows closely the method employed by Itzykson and 
Zuber to derive their well-known formula for the angular integral in hermitian 
matrix models. 

First we note, following that the partition function on the cylinder can easily 
be written in terms of the one on the disk, thanks to the invariance under area- 
preserving diffeomorphisms. One just has to deform the circle into a rectangle, 
decompose the Wilson loop in the product of four Wilson lines, and then identify 
two opposite sides. Consider then the SU{N) kernel on the disk, as given by |jTO[ 
before periodicity is imposed: 



/Ci(0.;t)=Ar(t)nf^^exp 



i<j 2 sin 2 
We want to calculate the integral 



2t 4" 



(45) 



lC2{^\,(Pi-t)= f dh]CiiAihA2h-^;t), (46) 

JSU{N) 

where Ai = diag{e^'^'^), A2 = diag{e'^'^'^), and /Ci is a function of the eigenvalues e'^'^u 
of the unitary matrix U = AihA2h~^. To this end, consider an arbitrary solution of 
the heat equation on the group, subject to the condition that it should be a class 
function, and a symmetric function of the eigenvalues e^'^u. Denote this function by 
/(0i;;t). Then 

f{<Pu;t) = I dglC,{Ug-\t)f{g,Q) 

JSU(N) 

dglC{U,g;t)fig;0). (47) 

SU{N) 

Following standard techniques, we integrate first over the eigenvalues of g, obtaining 
/(0i;; t) = CN I dS I dA,j\Ag)K{Au, SAgS-'; t)/(A,; 0), (48) 

JSU{N) J 

where 

J(A,)=.n2sin^^^^^. (49) 
i<j ^ 

Consider now the antisymmetric function of Ay 

^{<Plj,t)^A{Au)fi(P\;,t). (50) 
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We find 



= / dKgi[ct>%Q)K{Kg,Ku-.t), (51) 
where the kernel for ^ is 

/C(Ag, At/; t) = CA(A,) A(Ac;) / dSlC{Ku. SAgS-'; t). (52) 

JSU(N) 



Since f{4>if,t) is a solution of the heat equation , ^{(f)\j,t) obeys 

^ = Ai^u)^ufi^u;t), 



(53) 



where Ajj is the "radial" part of the Laplace-Beltrami operator on the group man- 
ifold, given by [|10| 

(54) 

with Rn = j^N^N"^ — 1). The function C,{4>jj,t) thus obeys the simple diffusion 



equation 



(55) 



For i?7v = 0, ^ is a solution of a diffusion equation which is totally antisymmetric in 
its arguments (pi, so the corresponding kernel must be the antisymmetric gaussian 



^(01,02;^) 



W^E(-l)''exp 
p 



(27rt)^ AT! 



_ ,(2) n2 



i2Trty 



det 



exp I 



Equating eq. ( |56D and eq. (|5^) we find that our integral is given by 



(56) 



SU{N) 



dh}C{Ai,hA2h-';t) 



JiA\)JiAl 



■det 



exp 



2t 



(57) 



and the only modification due to the R^ term is an extra factor of e^^*/^, which is 



also present in |10 
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Note 

While writing this paper we received two independent papers |]15|, |T^ with 



new results on QCD2, that partially overlap our work. In particular the work of 
Minahan and Polychronakos |T5[ was of some help to us in the interpretation of 



the grand canonical partition functions. 
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